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Abstract 

We present new sets of fragmentation functions for charged pions, charged kaons, 
and protons, both at the leading and next-to-leading orders. They are fitted to the 
scaled- momentum distributions of these hadrons measured in e + e" annihilation on 
the Z-boson resonance at CERN LEP1 and SLAC SLC. These data partly come 
as light-, charm-, bottom-quark-enriched and gluon-jet samples, which allows us to 
treat all partons independently, after imposing the SU(2) flavour symmetry rela- 
tions. In order to gain sensitivity to the scaling violation in fragmentation, we also 
include data from SLAC PEP, with center-of-mass energy ^/s = 29 GeV, in our fits. 
This allows us to also determine the strong-coupling constant, with a competitive 
error. LEP1 data on the longitudinal cross section as well as DESY DORIS and 
PETRA data at lower energies nicely agree with theoretical predictions based on 
our fragmentation functions. 

PACS numbers: 13.65.+i, 13.85.Ni, 13.87.Fh 



1 Introduction 



The inclusive production of a single charged hadron, h, in the annihilation process 

e+e" — (7, Z) -> h + X, (1) 

has been measured at many different e + e~ colliders over a wide range of center-of-mass 
(CM) energies, y/s, between 3 and 183 GeV Here, h may either refer to a specific 

hadron species, e.g., pion, kaon, or proton, or to the sum over all charged-hadron species. 
A large amount of precise data has now become available from various experimental 
collaborations at CERN LEP1 and SLAC SLC, who started taking data several years 
ago. The process (D) is particularly suitable in order to study the fragmentation of quarks 
and gluons into hadrons. The information contained in fixed-target, hadron-collider, and 
ep-scattering data is less useful, since it is obscured by theoretical uncertainties from the 
parton distribution functions and the choice of factorization scales connected with the 
initial states. 

The partonic cross sections pertinent to process (|l|) can be entirely calculated in per- 
turbative QCD with no additional input, except for the strong-coupling constant a s . They 
are known at next-to-leading order (NLO) f3[ and even at next-to-next-to-leading order 
(NNLO) |4j]. The subsequent transition of the partons into hadrons takes place at an 
energy scale of the order of 1 GeV and can, therefore, not be treated in perturbation 
theory. Instead, the hadronization of the partons is described by fragmentation functions 
(FFs), D^(x,Q 2 ). Their values correspond to the probability that the parton a, which 
is produced at short distance, of order 1/Q, fragments into the hadron h carrying the 
fraction x of the momentum of a. In the case of process (0), Q 2 is typically of order s. 
Given their x dependence at some scale Qq, the evolution of the FFs with Q 2 may be 
computed perturbatively. At present, the evolution equations are only known through 
NLO. Consistency requires that we do not include the NNLO corrections to the partonic 
cross sections [SJ in our analysis of process ([T]) via unpolarized photons and Z bosons. 

The theoretical analysis of the fragmentation process in e + e~ annihilation is of interest 
for several reasons. First, it allows us to test QCD quantitatively within one experiment 
observing hadrons at different values of \/s and to determine the strong-coupling constant, 
a s (Q), to be compared with determinations from other observables and/or processes. Sec- 
ond, since according to the factorization theorem the FFs are independent of the process 
in which they have been determined, they can be used for quantitative predictions of in- 
clusive single hadron cross sections in other processes, like pp, ep, jp, and 77 scattering. 
Third, the exact knowledge of the FFs, in particular for identified hadrons, may help to 
elucidate the fundamental features of hadronization and to constrain models used for the 
calculation of complete final states in various high-energy reactions. 

After the pioneering leading-order (LO) analyses of pion, kaon and charmed- 

meson FFs in the late 70s, there had long been no progress on the theoretical side of 
this field. In the mid 90 's, NLO FF sets for 7r° [Q, tt^, K , and rj mesons || have been 
constructed through fits to data of e + e~ annihilation, mostly generated with Monte Carlo 
event generators. 
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In 1994/95, two of us, together with Binnewies, determined the FFs of quarks and 
gluons into and mesons at LO and NLO from genuine experimental data, adopting 
two different strategies |lQ,|H[. In the first analysis flOfl , we included in the fit only the data 
on 7T and production taken by the TPC/Two-Gamma Collaboration |12] at SLAC 
PEP, with energy = 29 GeV. These data combine small statistical errors with fine 
binning in x and are thus more constraining than data collected by other experiments in 
the energy range from 5.2 GeV to 44 GeV. Charged-pion and -kaon data at other energies 
were only used for cross checks. Furthermore, theoretical predictions obtained using these 
FFs were compared with data on the inclusive production of unidentified charged hadrons 
taken at the PEP, DESY PETRA, and KEK TRISTAN colliders making an assumption 
on the p/p contribution based on information extracted from the TPC data |12[]. In 1994, 



when the first analysis |T0[ was performed, data with identified quark flavours did not 
exist. Therefore, we had to make the assumption that the s, c, and b (d, c, and b) quarks 
fragment into ^ {K^) mesons in the same way. 

In our second analysis we could remove this assumption because new data on 
charged-hadron production with partial flavour separation had been released by the 
ALEPH Collaboration [13| at LEP1 (see also Refs. @,|14|,[15| ) . Apart from the usual 
data sample without flavour separation, they also presented light- and 6-quark-enriched 
samples. Here and in the following, we consider u, d, and s quarks as light. In addition, 
we included in our fits accurate ir^ and data from ALEPH jl6| and information on 
tagged three-jet events from OPAL [|T7J at LEP1, which constrained the gluon FFs. In 
order to describe the production of charged hadrons in terms of the ir^ and FFs, 
we needed some information on the fragmentation of quarks and gluons into protons and 
antiprotons. To that end, we introduced a function f(x) which parameterizes the ratio of 
the p/p and 7r production cross sections as measured by ALEPH ||TB| . 

It is clear that this analysis JTTJ] suffered from the lack of specific data on the frag- 
mentation of tagged quarks and gluons to n^, K ± , and p/p hadrons, rather than just to 
charged hadrons of unidentified species. This lack has been remedied in 1998 by the ad- 
vent of new data from DELPHI pf at LEP1 and SLD @ at SLC. Besides other things, 
SLD [^] measured the differential production cross section, as a function of the scaled 
momentum x = 2p/^fs, for the identified hadron species n^, K^, and p/p in Z-boson 
decays to light, c, and b quarks and in such decays without flavour separation. Similarly, 
DELPHI [|18| presented n^, K , and p/p production data, covering the full momentum 
range up to 45.6 GeV, discriminating between Z — > uu, dd, ss, Z — > bb, and Z — > qq 
events. There is also new specific information that allows us to better control the gluon 
FF. In fact, new gluon-tagged three-jet data of inclusive charged-hadron production was 
released by ALEPH @ and OPAL In the OPAL |T| analysis, the gluon jets were 
identified using the method proposed by Gary |22]| . This method is based on rare events of 
the type e + e~ — > qqg inc \ in which the q- and g-quark jets appear in the same hemisphere of 
a multihadronic e + e~ annihilation event. The object gi nc i, which is taken to be the gluon 
jet, is defined by all hadrons observed in the hemisphere opposite to that containing the 
q- and g-quark jets. This method allows for the extraction of the gluon FF with as little 
theoretical bias as possible. It is supposed to be superior to the earlier OPAL [17] mea- 
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surement of the gluon FF based on a jet-finding algorithm to define the gluon jets, which 
entered the BKK analysis [[n]. The new OPAL ]21| data set comprises more points than 



the previous one fPTfl , and the energy of the gluon jets is now well defined, with a mean 
value of E 



'jet 



40.1 GeV. 



It is the purpose of this work to make full use of all the recent data from ALEPH 



2g, DELPHI Jl8|, OPAL |Hf, and SLD |9[ in order to construct new sets of LO and 



NLO FFs without imposing relations between the FFs of different flavours, except that 
we identify the FFs of u and d quarks into ir^ mesons and those of u and s quarks into 
mesons, as in Refs. [ |10|JTT|] . In the case of protons and antiprotons, we fix the <i-quark 
FFs to be half of the w-quark FFs. For completeness, we also include in our fits the older 
tt^, K , and p/p production data without flavour separation from ALEPH 
also entered the BKK analysis 



16 , which 



TTJ. In order to be sensitive to the running of a s 



i.e., 



to be able to fit the asymptotic scale parameter 
12 1, as in Refs 



A^, we also include the TPC data 



There is also precise data on charged-hadron production with 



flavour separation from ALEPH jjTJ|,[n|, who provided light-, o, and 6-quark-enriched 
samples, and from OPAL | 23| , who provided light- and 6-quark-enriched samples. To be 
on the conservative side, we exclude this data from our fits because the contribution due 
to charged hadrons other than 7r ± , K , and p/p is unknown and may be comparable 
or larger than the experimental errors. In fact, it turns out that this data yields rather 
sizeable y 2 values relative to the theoretical predictions obtained from our FFs. We also 
test our FFs against data from e + e~ colliders with lower CM energies, namely, from DESY 
DORIS [24,25] and PETRA [26,27]. Furthermore, we study the longitudinal cross section 
for the production of charged hadrons on the Z-boson resonance and compare the results 
with data from ALEPH 



DELPHI 28 



and OPAL 

This work is organized as follows. In Section 2, we summarize the data included in 
our fits, outline the theoretical framework, describe our fitting procedure, present our 
FFs, and discuss their goodness. Furthermore, we check our FFs against e + e~ data which 
are not included in our fits. In Section 3, we investigate for charged-hadron production 
how well our FFs satisfy the momentum sum rule, which we did not impose on our fits. 
Furthermore, we examine to what extent our FFs also account for the longitudinal cross 
section of charged-hadron production on the Z-boson resonance, and how our NLO gluon 
FF compares with the one of the BKK set [ 1 1] and with the one recently measured by 



DELPHI |[30||. Our conclusions are summarized in Section 4. 



2 Analysis and Results 



We start by summarizing the experimental data which are used for our fits or for com- 
parisons. We include in our fits the 7r ± , K , and p/p data with flavour separation from 
DELPHI |] 
PHI 



and SLD [[0| and those without flavour separation from ALEPH |I6| , DEL- 

Furthermore, we use the charged-hadron data with 



T8J, SLD 0, and TPC |2 



flavour separation from DELPHI and those without flavour separation from SLD |TJ| . 
In order to constrain the gluon FFs, we employ the gluon-tagged three-jet data without 



4 



Table 1: CM energies, types of data, and Xdf values obtained at LO and NLO for various 
data samples. Samples not used in the fits are marked by asterisks. 



[GeV] 


Data type 


in NLO 


Xb F in LO 


29.0 


(T* (all) 
a K (all) 
a p (all) 


0.64 [ 
1.86 [ 
0.79 [ 


12 
12 
L2 




0.71 [ 
1.40 | 
0.70 [ 


12 
12 
12 




91.2 


a h (all) 

a' 1 (uds) 
a h (c) 
a h (b) 


1.28 [ 
1.32 [ 
0.20 [ 

0.43 [ 


18 
19 
18 

18 


48.1 |15(1* 14.4 

89.4 @* 5.10 [ 
80.1 111* 0.51 [ 
221 (lj* 6.63 


23f 

n* 

2l* 
23J* 


1.40 [ 
1.44 [ 
0.20 | 

0.41 


18 
19 
18 

E 


53.1 |15|* 15.8 [ 

92.1 @* 4.60 [ 
58.9 JTJ* 0.45 
220 p4|* 6.33 [ 


23 



23 
23 
23 


1* 

1* 
1* 
1* 




ex 77 (all) 
a* (uds) 
a' (c) 

CT- (b) 


1.28 [ 


16 


0.58 [ 
0.72 

0.57 | 


18 
18 

18 


] 3.09 [ 
] 1.87 | 
1.36 [ 
] 1.00 


1GJ 

m 
lg 
il 


1.65 


16 


0.60 [ 
0.73 | 

0.58 [ 


18 
18 

18 


3.13 [ 
2.17 [ 
1.16 [ 
0.99 [ 


19 
19 
19 

i9 






a K (all) 
<x K (uds) 
a K (c) 
<7* (b) 


0.30 |6| 0.86 [ 
0.53 

0.14 | 


18 
18 

18 


] 0.44 [ 
j 0.65 | 
2.11 [ 
1 1-21 


1ST 
m 
lg 
il 


0.32 |16|] 0.79 [ 
0.60 | 

0.14 [ 


18 
18 

18 


0.45 [ 
0.64 | 
1.90 | 
1 1-23 | 


19 
19 
19 
19 


1 




erf (all) 
cr p (uds) 
^ (c) 
a* (b) 


0.93 |6| 0.09 [ 
0.11 

0.56 | 


18 
18 

EL8 


] 0.79 [ 
] 1.29 [ 
0.92 [ 
] 0.97 | 


lpj 

13 
il 
il 


0.80 |lj 0.06 [ 
0.14 | 

0.62 [ 


18 
18 

EL8 


0.70 [ 
1.28 [ 
0.89 [ 
0.89 | 


19 
19 

is 

19 




E iet [GeV] 








26.2 
40.1 




1.19 [E0 
1.03 [gl 




1.18 pi 
0.90 || 



hadron identification from ALEPH |2(J and OPAL | 21| . The charged-hadron data with- 
out and with flavour separation from ALEPH [|T4| , [T5| and OPAL |H| are only used for 
comparison. All these data sets are summarized in Table 1. 

The key observable of the experimental analyses and our study is the scaled-momentum 
distribution normalized to the total hadronic cross section, (l/a tot )dcr h /dx. In general, 
the scaled momentum is defined as x = 2p/^, where p is the three-momentum of the 
observed hadron in the CM frame. In three-jet events, one often defines x = p/Ej Ct , where 
i?j et is the gluon-jet energy in the CM frame. By charge-conjugation invariance, the e + e~ 
cross sections for ir + , K + , and p production should be the same as those for 7r~, K~ , and 
p production, respectively. Therefore, one usually sums over both charges, which is also 
true for the production of unidentified charged hadrons. In turn, also our FFs refer to 
the sums of particles and antiparticles. 

The LO and NLO formalisms for extracting FFs from e + e~ data were comprehensively 
described in Refs. ||T()| , |TTf and will, therefore, not be reviewed here. The NLO formula for 
da h /dx was originally obtained in Ref. and is given in Eq. (3) of Ref. Deviating 
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from Refs. [|K],|ri|], we calculate <7 to t including QCD corrections up to 0(a 2 s ). We work 
in the MS renormalization and factorization scheme and choose the renormalization scale 
ft and the factorization scale Mf to be fj, — Mf = a/s, except for gluon-tagged three- 
jet events, where we put \i = Mf = 2E- ]Ct . For the actual fitting procedure, we use x 
bins in the interval 0.1 < x < 1 and integrate the theoretical functions over the bins 
width as is done in the experimental analyses. The restriction at small x is introduced 
to exclude events in the nonperturbative region, where mass effects and nonperturbative 
intrinsic-transverse-momentum (fcy) effects are important and the underlying formalism 
is insufficient. We parameterize the x dependence of the FFs at the starting scale Qo as 

D h a (x,Ql) = Nx a (l-xf. (2) 

As in Ref. |TT|, we use Q = \/2 GeV for a = u,d,s,g, Qo = m(r] c ) = 2.9788 GeV for 
a = c, and Qo = m 00 = 9.46037 GeV for a = b. We impose the conditions 



DZ (x, Qi) = Dl (x, Q 
Df{x,Ql) = Df{x,Ql), 

D^(x,Ql) = 2Df(x,Ql), (3) 

which are suggested by the valence quark composition of the respective hadrons. Equa- 
tions (H) are preserved by the Q 2 evolution. The constraint on the K FFs should be 
slightly violated by the mass difference between the u and s quarks ||. In fact, the 
s — > K~ transition should happen more frequently than the u — > K~ one because less en- 
ergy is needed for the creation of a uu pair from the vacuum than for a ss pair. However, 
we shall neglect this subtlety in our analysis because the presently available experimental 
information is not rich enough to resolve such minor differences. For all the other FFs, 
we treat N, a, and 8 as independent fit parameters. In addition, we take the asymptotic 
scale parameter appropriate for five quark flavours, as a free parameter. Thus, we 
have a total of 46 independent fit parameters. The quality of the fit is measured in terms 
of the x 2 value per degree of freedom, Xdfi f° r a ^ selected data points. Using a multi- 



dimensional minimization algorithm [5T|, we search this 46-dimensional parameter space 
for the point at which the deviation of the theoretical prediction from the data becomes 
minimal. 

The Xbf va hies achieved for the different data sets used in our LO and NLO fits may 
be seen from Table 1. Most of the Xbf values lie around 1 or below, indicating that the 
fitted FFs describe all data sets within their respective errors. In general, the Xbf values 



come out in favour of the DELPHI [T8|] data. The overall goodness of the NLO (LO) fit is 



given by Xbf =0.98 (0.97). In Table 1, we have also specified the Xdf values for those data 



sets which we have discussed above, but not included in the fits [14,15,23] . Some of these 
sets come out with much higher Xdf values than those included in the fits, which will be 
discussed in detail later. The values of the parameters N, a, and 8 in Eq. (|) resulting 
from the LO and NLO fits are collected in Table 2. They refer to the fragmentation of the 
partons a — u,d, s, c, b, g into the hadrons h = 7r ± , K ± ,p/p, where the sum over particles 
and antiparticles is implied. 
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Table 2: Values of N, a, and (3 in Eq. (fj) resulting from the NLO fit. The numbers given 
in parentheses refer to the LO fit. 



Hadron 


Flavour 


AT 

iv 




Q' 




p 


7T 


u = d 


H A A Q ( C\ t^A(Z\ 

U.44o {v.owj 


1 A Q 


I 1 A 7\ 


n mo 


(1.02) 




s 


1D.0 [ZZ.o) 


ft 1 QQ 


(U.LZ ( ) 


K on 


(R 1 /I "\ 




c 


b.17 (8.76) 


—0.53b 


(—0.386) 


5.60 


(5.62) 







n oc;n (r\ qi i \ 
U.zoy (U.oii J 


i no 

— i.yy 


( 1 OQ\ 

\ — l.yo j 


O.OO 


[oAl ) 




n 

y 


3 73 (6 05) 


-0.742 


(-0.714) 


2.33 


(2.92) 




u = s 


0.178 (0.259) 


-0.537 


(-0.619) 


0.759 


(0.859) 




d 


4.96 (5.38) 


0.0556 ( 


-0.00321) 


2.80 


(3.08) 




c 


4.26 (5.18) 


-0.241 


(-0.178) 


4.21 


(4.30) 




b 


1.32 (1.57) 


-0.884 


(-0.841) 


6.15 


(6.01) 




9 


0.231 (0.0286) 


-1.36 


(-2.94) 


1.80 


(2.73) 


P/P 


u = 2d 


1.26 (0.402) 


0.0712 


(-0.860) 


4.13 


(2.80) 




s 


4.01 (4.08) 


0.173 ( 


-0.0974) 


5.21 


(4.99) 




c 


0.0825 (0.111) 


-1.61 


(-1.54) 


2.01 


(2.21) 




b 


24.3 (40.1) 


0.579 


(0.742) 


12.1 


(12.4) 
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1.56 (0.740) 


0.0157 


(-0.770) 


3.58 


(7.69) 



Since we included in our fits high-quality data from two very different energies, namely 
29 and 91.2 GeV, we are sensitive to the running of a s (n) and are, therefore, able to 
extract values of A||. We obtain A|| = 213+^ MeV at NLO and A|| = 88+^ MeV at 
LO. This corresponds to a s (M z ) = 0.1170+^ at NLO and a s (M z ) = 0.1181+^ at 

LO. Here, the errors are determined by varying A^- in such a way that the total Xdf * s 
increased by one unit if all the other fit parameters are kept fixed. They do not include 
theoretical errors, which may be estimated through variations of the overall scale /i = Mf. 
Our NLO values for A^ and a s (Mz) perfectly agree with those presently quoted by the 
Particle Data Group as world averages, A^ = 212+^ MeV and a s (M z ) = 0.1185±0.002, 

respectively | 32|| . We note, however, that in Ref. [32] the value of A^ is obtained from 
a s (Mz) through the next-to-next-to- leading-order formula for a s (fi), while we use the LO 
and NLO formulas, depending on the order of our analysis, as is required by consistency. 
Since the errors on a s (jj) and A^ quoted in Ref. f32| result from an average of different 
kinds of measurements, they are considerably smaller than those obtained here from a 
single type of experiment. In this context, we remark that, throughout our analysis, we 
assume nonperturbative power corrections proportional to 1/Q to be negligible in the 
energy range relevant for our fits. 

The goodness of our LO and NLO fits to the ALEPH 1TJ,|20[, DELPHI (T§, OPAL 



2j| , and SLD Jl9| data may be judged from Figs. 1-5. In Figs. 1-4, we study the cross 
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section (l/a to t)da h /dx of charged-hadron, 7r ± , K ± , and p/p production, respectively, at 
yfs = 91.2 GeV as a function of x. In Fig. 1, the contributions from all quark flavours 
and the gluon are included, while Figs. 2-4 refer to the fragmentation of light, c, and b 
quarks, respectively. As in the experimental analyses |T4] , p!5| , p!8| , p!9| , p3| , we have to arrange 
for the latter three contributions to add up to the first one. Therefore, we distribute the 
gluon contribution among the various quark contributions according to their electroweak 
coupling strengths to the neutral current. In this way, the gluon radiation off a quark 
line of a given flavour is accounted towards the contribution of this quark flavour. The 
theoretical results are compared with the respective data from ALEPH |16| in Fig. 1, with 
that from DELPHI 1181 in Figs. 1, 2, and 4, and with that from SLD pf in Figs. 1-4. 



fL8| in Figs. 1, 2, and 4, and with that from SLD 
We observe that, in all cases, the various data are mutually consistent with each other 
and are nicely described by the LO and NLO fits, which is also reflected in the relatively 
small Xdf values given in Table 1. The LO and NLO fits are almost indistinguishable in 
those regions of x, where the data have small errors. For large x, where the statistical 
errors are large, the LO and NLO results sometimes moderately deviate from each other. 

In Fig. 5, we compare the ALEPH |2l| and OPAL |2IJ] measurements of the gluon FF 
in charged-hadron production, which comprises five and twelve data points, respectively, 
with our LO and NLO fit results. The scale choice is fi = Mf 



2E jeU with Ejet 



26.2 



|2(J and 40.1 GeV pi] , respectively. The data are nicely fitted, with Xt>f values of order 



unity, as may bee seen from Table 1. By the same token, this implies that the data are 
mutually consistent. 

As may be seen from Table 1, the ALEPH [|l4|,[il| and OPAL |23| data on charged- 
hadron production with and without flavour separation, which are excluded from our fits, 
lead to rather sizeable Xdf values relative to the theoretical predictions based on our FFs. 
They are obviously inconsistent with the DELPHI 1TS[ and SLD [fill data included in 
our fits. In particular, the flavour-enriched charged-hadron samples from ALEPH 



give Xdf values of order 100. The comparison of these ALEPH fll^JIBI and OPAL [23] 
data with our LO and NLO predictions is also illustrated in Fig. 6. We should point out 
that the ALEPH [T^JT^,[16| data possess the following property. The sum of the three 
flavour-tagged cross sections of charged-hadron production JTJI, which agrees with the 
charged-hadron cross section without flavour separation quoted in Ref. ||15|| , appreciably 
overshoots the sum of the n^, K ± , and p/p cross sections without flavour separation, 
which stem from an independent measurement [|16 



The latter three data sets [16 



are 

included in our fits and lead to quite acceptable Xdf values, as may be seen from Table 1. 
We remark that the charged-hadron data with flavour separation from ALEPH [14] served 
as one of the main inputs for the BKK analysis [111 ]. We attribute this difference to the 
contribution from charged hadrons other than ir ± , K , and p/p hadrons, which may be 
comparable or larger than the relatively small errors on the charged-hadron data p~4 JT5 1 . 
We assume that this interpretation is also valid for the OPAL data |]2"3" |. On the contrary, 
the DELPHI 



and SLD [19 



analyses are based on the assumption that the sum of 
the it^, K ± , and p/p data exhaust the full charged-hadron data. In want of a separate 
measurement of the excess in cross section due to high-mass charged hadrons, we decided 
to proceed as in Refs. [|i"8|.|T9"|. 
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Table 3: CM energies, types of data, and Xdf values obtained at LO and NLO for various 
pre-LEP/SLC data samples not included in our fits. 



[GeV] 


Data type 


xl F in NLO 


Xr> F in L0 


5.4 


(all) 
a K (all) 


3.10 [24 
1.80 [24 


1 
1 


3.02 [ 
2.33 [ 


24] 
24] 


9.98 


a w (all) 
a K (all) 


3.27 [25 
3.21 [25 


1 
1 


2.78 [ 
2.81 [ 


25] 
25] 


22.0 


aP (all) 


1.29 || 


1.50 || 


34.0 


a* (all) 
a K (all) 
a p (all) 


0.80 [27 
0.31 [27 
0.58 [27 


1 
1 
1 


0.88 [ 
0.37 [ 
0.46 [ 


27] 
27] 
27] 



We now report on comparisons with pre-LEP/SLC data on 7r ± , K , and pp production 
which we did not include in our fits. Specifically, we consider the 7r ± and data from 
DASP |2| (y/s = 5.2 GeV) and ARGUS || (y/s = 9.98 GeV) at DESY DORIS and from 
TASSO H (V£ = 34 GeV) at PETRA, and the pjp data taken by TASSO at ^ = 22 
|26| and 34 GeV |[27|| . There is no separation into flavour-enriched or three-jet samples 
for this data. The resulting Xdf values are summarized in Table 3. As in the fits, they 
are evaluated only taking into account the data points with x > 0.1. The comparison 
of the 7r ± , K ± , and pjp data with our LO and NLO predictions is visualized in Figs. 7- 
9, respectively. For reference, also the corresponding TPC [|12| and SLC [I9|] data are 
included in these comparisons. From Table 3 and Figs. 7-9, we learn that the DASP f24| , 
ARGUS 



25 , and TASSO 



data agree quite well with our LO and NLO predictions. 



This concludes the description of our fitting procedure and of the quality of the fits. 



3 Applications 

We now investigate if our FFs satisfy the momentum sum rules. Since a given parton a 
fragments with 100% likelihood into some hadron h and momentum is conserved during 
the fragmentation process, we have 

1 dxxD h a (x,Q 2 ) = l, (4) 



o 

for any value of Q 2 . As is well established experimentally, even for the flavour-enriched 
data samples [|1!|,|]1|, the sum of the n ± , K ± , and pjp production cross sections agrees, to 
good approximation, with the one of charged-hadron production. Thus, insertion of our 
ir^, K^, and p/p FFs in Eq. (Q) exhausts the charged-hadron contribution. As for the 
neutral-hadron FFs, which also enter Eq. (H), we make the following assumptions which 
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are suggested by SU(2) flavour symmetry. As in Ref. JTT[] , we approximate 



Furthermore, we assume that 



i,s,c,b,g\ 

D n J n (x,Q 2 ) = ±D^(x,Q 2 ) 

D n d /n (x,Q 2 ) = 2Df(x,Q 2 ) 
l,n (r O 2 ) - D p/P 

:,c,b,g\ x ' V J — ^ s,c,b,g 



D:!l Q (x,Q 2 ) = D^ b (x,Q*). (6) 



Here, Df(x,Q 2 ), Df/ R °(x,Q 2 ), etc. refer to the sums of the FFs for the individual 
hadron species. Equations (|5|) are supported by e + e~ data of pion pLB3[ and kaon 
production, except for very small x. To be on the conservative side, we take the lower 
limit of integration in Eq. (f|) to be 0.05 rather than 0, so that the FF parameterizations 
are not used too far outside the x range which was selected for the fits. This particular 
choice may be motivated by observing that the n ± , K ± , and p/p data sets in Figs. 7-9 
are well described by the theoretical predictions down to x values of order 0.05. Since we 
leave out the x range below 0.05, our results for the left-hand side of Eq. (P are expected 
to be somewhat below unity. A more sophisticated approach would be to improve 
the description of the low-x region by modifying our ansatz for the FFs according to the 
so-called modified leading logarithmic approximation (MLLA) |p6| . 

In Table |], we list the results obtained with our LO and NLO FFs for Q = y/2, 4, 
10, 91, and 200 GeV. The accuracy of our results is limited by the uncertainties in our 
assumptions (|5]) and (|6|), in particular for the lower values of Q. In the case of the d quark, 
we find values slightly larger than unity at small values of Q, both in LO and NLO. On the 
other hand, the s-quark values are considerably smaller than unity, in particular in NLO. 
The gluon FF gives values appreciably larger than unity at small values of Q, both in LO 
and NLO. We believe that the imbalance between the d- and s-quark results is due to our 
limited knowledge of the relative magnitude of the valence- and light-sea-quark contents 
of the 7r ± , K^, and p/p hadrons. Since the u, d, and s quarks are always combined in 
the available data, there are no individual constraints on their FFs. In order to trace the 
source of this imbalance, let us consider the ir^ 1 and K ± contributions to the left-hand 
side of Eq. (HI) at Q = \p2 GeV for a = u,d, s. In the case of 7r ± , we have 



f dxxDZ%x,Ql) = UAl, 

J0.05 



dxxDl (x,Q5) = 0.23, (7) 

0.05 

which nicely exhibits the valence-like character of the u and d quarks in contrast to the 
sea-like character of the s quark. In the case of K^, however, we find 



f l dxxD«f(x,Q 2 ) = 0.19, 

J0.05 
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Table 4: Left-hand side of Eq. at NLO for a = u,d,s,c,b,g and Q = y/2, 4, 10, 91, 
200 GeV. We sum over h = n^, it , K ± , K°, K°,p,p, n, n and integrate over 0.05 < x < 1. 
The numbers given in parentheses are evaluated with our LO set. 



a 


Q [GeV] 




V2 


4 


10 


91 


200 


u 


0.96 (1.13) 


0.96 (1.06) 


0.92 (1.00) 


0.82 (0.87) 


0.78 (0.83) 


d 


1.05 (1.14) 


1.05 (1.07) 


1.00 (1.01) 


0.89 (0.88) 


0.85 (0.84) 


s 


0.60 (0.82) 


0.68 (0.82) 


0.67 (0.78) 


0.62 (0.69) 


0.60 (0.67) 


c 




0.91 (0.97) 


0.88 (0.93) 


0.78 (0.81) 


0.75 (0.78) 


b 






0.73 (0.80) 


0.65 (0.69) 


0.62 (0.66) 
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1.33 (1.85) 


1.15 (1.32) 


1.00 (1.10) 


0.73 (0.74) 


0.67 (0.67) 



0.05 



doc xD ^ 



(x,Q 



0.25, 



which shows that the d quark does not behave sea-like, contrary to expectations. We have 
thus identified the non-sea-like behaviour of the d quark in the K ± mesons as the source 
of the violation of the sum rule in Table 4. 

The analysis of the the longitudinal cross section offers a unique opportunity to test 
the gluon FF. The total cross section of inclusive hadron (h) production, which we have 
considered so far, may be decomposed into a transverse (T) and a longitudinal (L) compo- 
nent, da h /dx = da\ / dx + da\ / dx , as indicated in Eq. (3) of Ref. [|TTJ. As is well known, at 
0(aP s ), only transversely polarized photons and Z bosons contribute to the cross section. 
The longitudinal cross section da\jdx first appears at 0(a s ). Thus, one needs to know 
it through 0(a 2 s ) in order to test the gluon FF in NLO. The 0(a 2 s ) correction to da^/dx 
has been calculated in Refs. In the case of inclusive hadron production, da^/dx is 

given by a convolution of the longitudinal partonic cross sections with the corresponding 



FFs and has the following form |p7| : 



1 do\ 



o"tot dx 



mi"- 



JSf V-z / °tot dz y \z 1 ) cr to t dz 



i=i kz 



cr tot dz 



Nf 

i=i 



1 da[ tg 



dtot dz 



(9) 



Here, Nc = 3 and Np = 5 denote the numbers of colours and active quark flavours, 
respectively, and Q gi (s) and Qq.(s) represent the effective electroweak couplings of the 
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quarks to the photon and Z boson. The latter are defined as 



_ , , 2 2 Sis — Ml 

Q qi {s)^e e e +2e e v e e qi v t 



s - Ml) 2 + Mfr| 



:s-MZ) 2 + M*TV 



l z L z 



Np 



s - M 2 .) 2 + M' 2 7 T 



2r2 



(10) 



z- 1 - z 



where Mz and are the mass and width of the Z boson, respectively, ef is the electric 
charge of the fermion / in units of the positron charge, Vf — (t 3 j — 2ej sin 2 Oy^j / (2 sin^ 
cos9 w ) and aj = T 3 j / (2 sin 8 W cos 9 W ) are its vector and axial- vector couplings to the Z 
boson, respectively, T^j is the third component of isospin of its left-handed component, 
and 8 W is the electroweak mixing angle. In the limit s <C Mf, only the QED contribution 



survives, and we have Q qi 



HI 



and Q q _ — * 0. The labels £ and NS in Eq. (H) indicate 



the flavour singlet and flavour non-singlet contributions, which are related to the FF 
combinations 



(x, Mf) = £ [D h qi (x, M|) + D\ (x, M) 

M f) = K (*, M J) + D i (*, M f) - ( x ' M f) 



The contributions proportional to the colour factor CfTr = 4/8 are labeled F, which 
is to suggest that they vanish in QED due to Furry's theorem. They only arise from 
Z-boson exchange, provided one does not sum over all quark flavours belonging to one 
fermion generation. Furthermore, they are devoid of collinear singularities, so that mass 
factorization is not needed. We neglect these contributions, since they amount to less 
than 1% over the whole x range []37|. The partonic cross sections da^ q /dz, da^g/dz, and 



daf s q /dz through 0(a 2 s ) may be found in Refs. j§g7|. 

ALEPH 0, DELPHI f28|, and OPAL || have measured da h L /dx for charged-hadron 
production as a function of x without quark flavour separation. In addition, DELPHI 
p8| has also presented light- and 6-quark-enriched samples. In Fig. 10, these data are 
compared with the LO and NLO evaluations of Eq. ([|) with our FFs. As in our fits to 
data of da h /dx, we choose the scales to be equal to the CM energy yfs = 91.2 GeV, where 
the data have been taken. Although only data with x > 0.1 has been included in our fits, 
da^/dx is also well described down to x = 0.01. The Xdf values, evaluated from the data 
points with x > 0.1, are collected in Table 5. The Xdf va hies for the LO predictions are 
in average more than a factor of two larger than those for the NLO predictions. This is 



not surprising, since it was already found in Ref. |TT[ that the LO result falls short of the 
data by a factor of two. At that time, the common scale had to be reduced to 20 GeV 
in order to obtain agreement. We observe from Fig. 10 and Table 5 that this is no longer 
necessary for the NLO results. Obviously, the 0(a 2 s ) corrections H,|37|| are sufficient to 
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Table 5: CM energies, types of data, and Xdf values obtained at LO and NLO for the 
longitudinal cross section. 



[GeV] 


Data type 


xl F in NLO 




91.2 


a h L (all) 
a\ (uds) 
(b) 


1.36 | 


1 1.74 |H 0.49 H 

7.98 m 

0.51 H 


11.0 [[[ 


| 1.27 |2S] 7.94 
1.05 |2S] 
0.80 |2S] 


2S] 



produce a K factor that brings our predictions into good agreement with the data. Good 
agreement between the NLO prediction of da\jdx and the ALEPH |14| and OPAL [2E 



data was already found by Rijken and van Neerven [|J, who employed the BKK gluon FF 

The comparisons presented in Fig. 10 provide us with a very useful check of the gluon 
FF in the low- a; range, where the data have small errors. In particular, the good agreement 
in the case of flavour separation nicely confirms the relative importance of the light- and 
6-quark FFs relative to the gluon FF. Since the Xdf values obtained at NLO for the 
longitudinal cross section are comparable to those achieved in our NLO fit, we do not 
expect that the inclusion of the information on the longitudinal cross section in our NLO 
fit would significantly change the outcome. 

As is well known, the gluon FF only enters the prediction for da h /dx at 0(a s ), while 
at 0{oP s ) it only contributes indirectly via the Q 2 evolution. Therefore, we have used the 
ALEPH [2(J and OPAL [21] data on gluon-tagged three-jet events in order to constrain 
it. These data were already compared with the fitted gluon FF in Fig. 5. It is of interest 



to compare this gluon FF with the BKK |Tl one and the one recently extracted by 



DELPHI [[!(J from their three-jet events for scales in the range 6.5 < Mf < 28 GeV. 
In Ref. |30| , a parameterization, valid for this Mf range and for 0.15 < x < 0.75, was 
generated adopting an ansatz similar to the one introduced in Ref. [|TTJ. In Fig. 11, the x 
dependences of these three gluon FFs are studied for Mf = 10, 52.4, 80.2, and 200 GeV. 
The results for M f = 52.4 and 80.2 GeV are compared with the ALEPH and OPAL 



21] gluon-tagged three-jet data, which are included in our fits. From Fig. 11 the following 



observations can be made. Our new gluon FF is rather similar to the BKK one, especially 
for larger values of Mf. At Mf = 80.2 GeV, the DELPHI gluon FF is somewhat steeper 
than the new and BKK ones, except for large x values, and it agrees slightly better with 
the OPAL |2TJ data. On the other hand, our new gluon FF agrees best with the ALEPH 
|3 data at M f = 52 .4 GeV. At M f = 10 GeV, the DELPHI gluon FF has a flatter 
x dependence than the other ones and overshoots them, while at Mf = 200 GeV the 
situation is similar to the case of Mf = 80.2 GeV. 



13 



4 Conclusions 



The LEP1 and SLC experiments have provided us with a wealth of high-precision experi- 
mental information on how partons fragment into low- mass charged hadrons P,|T3*,I%, 15,1^,17 
The data partly come as light-, c-, and 5-quark-enriched samples without |p!4| , p!5| , p3| or with 
identified final-state hadrons (ti ± , K ± , and p/p) |18|JT9|| or as gluon-tagged three-jet sam- 
ples without hadron identification pO|j2T| , |30[] . This new situation motivated us to update, 
refine, and extend the BKK analysis |TT| by generating new LO and NLO sets of tt , 
K , and p/p FFs. The x distributions of the resulting FFs at their starting scales Qq 
are given by Eq. (|2|) with the parameters listed in Table 2. The Q 2 evolution is deter- 
mined by the timelike Altarelli-Parisi equations in the respective orders, LO or NLO, 
which are summarized in the Appendix of Ref. The evolution procedure in Mellin 

space is described in Ref. [f[l"1 . A FORTRAN subroutine which returns the values of 
the D^(x, Q 2 ) functions for given values of x and Q 2 may be downloaded from the URL 
http : / /www . desy . de/~poetter/kkp . html] or obtained upon request from the authors. 

By also including in our fits tt , K , and p/p data (without flavour separation) from 
PEP ||12|| ) with CM energy y/s = 29 GeV, we obtained a handle on the scaling violation 
in the fragmentation process, which allowed us to extract LO and NLO values of a s [Mz)- 
We found a s (M z ) = 0.1181 +J;JJ? at LO and a s (M z ) = 0.1170 +g;gg* at NLO, where the 
errors are experimental. These results are in perfect agreement with what the Particle 



Data Group currently quotes as the world average, a s (M z ) = 0.1185 ± 0.002 p2 



Our strategy was to only include in our fits LEP1 and SLC data with both flavour 

gluon-tagged three-jet samples with a fixed 



separation and hadron identification ||19| , |18 
gluon-jet energy 



20.21 



and the 7r ± , K^, and p/p data sets from the pre-LEPl/SLC 
era with the highest statistics and the finest binning in x ||12| . Other data served us for 
cross checks. In particular, we probed the scaling violation through comparisons with 
71^, K^, and p/p data from DORIS and PETRA, with CM energies between 5.4 and 
34 GeV p4] , p5| , p6| , p7| . Furthermore, we tested the gluon FF, which enters the unpolarized 
cross section only at NLO, by comparing our predictions for the longitudinal cross section, 
where it already enters at LO, with available data [|14],^|^|. Finally, we directly compared 
our gluon FF with the one recently measured by DELPHI in three-jet production with 
gluon identification as a function of x at various scales Q 2 . All these comparisons led 
to rather encouraging results. We also verified that our FFs satisfy reasonably well the 
momentum sum rules, which we did not impose as constraints on our fits. All these cross 
checks make us believe that our FFs should allow for a reliable description of inclusive 
charged-hadron production in all kinds of high-energy-collision experiments over wide 
ranges of x and Q 2 . The very-high-Q 2 regime will be accessible with the CERN LHC and 
a future e + e~ linear collider. Our FFs are bound to break down at very low values of x, 
where it becomes necessary to modify the ansatz (|j) according to the MLLA [36|. This 
is beyond the scope of the present analysis, but should eventually be investigated 



1 There is an obvious typographical error in the published version of Ref. p8[ , which was absent in 
the preprint version thereof. In the line before the last of Eq. (17), lnln(l — x) should be replaced by 
ln(l-x). 
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Note added 



After the submission of this paper, a preprint [3J| appeared which also presents NLO FFs 
for n ± ,K ± , and charged hadrons. These FFs are fitted to the ALEPH [0,0, SLD |TJ, 
and TPC [O] data assuming certain power laws in the light-quark sector. 
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Figure 1: Normalized differential cross section of inclusive hadron production at y/s = 
91.2 GeV as a function of x. The LO (dashed lines) and NLO (solid lines) fit results are 



compared with data from ALEPH [|TJ (triangles), DELPHI (TJ (circles), and SLD [|TJ 
(squares). The upmost, second, third, and lowest curves refer to charged hadrons, tt 1 * 1 , 
K ± , and p/p, respectively. Each pair of curves is rescaled relative to the nearest upper 
one by a factor of 1/5. 
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Figure 3: Same as in Fig. 1, but for the c-quark-enriched samples from SLD [|1J (squares). 
The last two points on the right belong to the charged-hadron sample. 
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Figure 5: Gluon FF for charged-hadron production as a function of x at Mf = 52.4 and 
80.2 GeV. The LO (dashed lines) and NLO (solid lines) predictions are compared with 

with Ejet = 26.2 GeV, (upper curves) and from OPAL 



three-jet data from ALEPH 
21 1, with £jet = 40.1 GeV (lower curves). The OPAL data and the pertinent predictions 



are rescaled by a factor of 1/100. 
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Figure 6: Normalized differential cross section of inclusive charged-hadron production at 
y'i = 91.2 GeV as a function of x. The LO (dashed lines) and NLO (solid lines) fit results 
are compared with data from ALEPH [|14] , [15[j (triangles) and OPAL |23[] (circles). The 
upmost, second, third, and lowest curves refer to the full, light-quark-enriched, c-quark- 
enriched, and 6-quark-enriched samples, respectively. Each pair of curves is rescaled 
relative to the nearest upper one by a factor of 1/5. 
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Figure 7: Normalized differential cross section of inclusive 7r ± production as a function of 
x at y/s = 5.2, 9.98, 29, 34, and 91.2 GeV. The LO (dashed lines) and NLO (solid lines) 



predictions are compared with data from DASP f2J, ARGUS |§, TPC @, TASSO |7|, 
and SLD |19|| . Upper curves correspond to lower energies. Each pair of curves is rescaled 
relative to the nearest upper one by a factor of 1/10. 
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Figure 9: Normalized differential cross section of inclusive p/p production as a function 
of x at y/s = 22, 29, 34, and 91.2 GeV. The LO (dashed lines) and NLO (solid lines) 
predictions are compared with data from TASSO (with \fs = 22 [^6| and 34 GeV [p7| ). 
TPC [|12| , and SLD |]19|| . Upper curves correspond to lower energies. Each pair of curves 
is rescaled relative to the nearest upper one by a factor of 1/10. 
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Figure 10: Normalized differential longitudinal cross section of inclusive charged-hadron 
production as a function of x at y/s = 91.2 GeV. The LO (dashed lines) and NLO (solid 
lines) predictions are compared with data from ALEPH @], OPAL and DELPHI @ 
without flavour separation and with light- and 6-quark-enriched samples from DELPHI 
|28| (in this order from top to bottom). Each pair of curves is rescaled relative to the 
nearest upper one by a factor of 1/30. 
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Figure 11: Gluon FFs for charged-hadron production as functions of x at Mf = 10, 
52.4, 80.2, and 200 GeV from the DELPHI J2TJ (dotted lines), BKK NLO fill (dashed 
lines), and new NLO (solid lines) sets. The predictions for Mf = 52.4 and 80.2 GeV 
are compared with the gluon-tagged three-jet data from ALEPH f2(J and OPAL |2l| , 
respectively. Upper curves correspond to lower energies. Each set of curves is rescaled 
relative to the nearest upper one by a factor of 1/50. 
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